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ABSTRACT 

CDM simulations predict that there are hundreds of lumps of with masses 
O i greater than lO^M© in the Milky Way halo. However, we know of only a dozen 

^ I dwarf satellites close to this mass. Are these lumps simply lacking in stars or is 

Q I there a fundamental flaw in our most popular cosmology? By studying the tidal 

| 7j I debris of known satellites we can potentially address this question. In this paper, 

Q\ I we quantify the the effects of the dark matter lumps on tidal tails. The lumps 

scatter stars in the tidal tails from their original orbits producing a distinctive 
>• I signature. We simulate debris evolution in smooth and lumpy halos potentials 

^ ' and use our simulations to motivate and test a statistical measure of the degree 

^ ' of scattering apparent in the angular position and radial velocity measurement 

^ I of debris stars — the "scattering index". We find that the scattering index 

^ ■ can in general distinguish between the levels of substructure predicted by CDM 

^ ' cosmologies and smooth Milky Way models, but that the sensitivity of the debris 

Ph' depends on the orientation of the parent satellite's orbit relative to the largest 

O ' lumps orbits. 

• We apply our results to the carbon star stream associated with the Sagittarius 

■ dwarf galaxy (Sgr) and find that these stars appear to be more scattered than we 

^ ^ ■ expect for debris orbiting in a smooth halo. However, the degree of scattering is 

^ ' entirely consistent with that expected due to the influence of the Large Magellanic 

c3 ■ Cloud, which is on an orbit that intersects Sgr's own. We conclude that the 

current data is unable to constrain CDM models. 

Nevertheless, our study suggests that future data sets of debris stars associ- 
ated with other Milky Way satellites could provide strong constraints on CDM 
models. 

Subject headings: Galaxy: halo — Galaxy: kinematics and dynamics — Galaxy: 
structure — cosmology: dark matter 
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1. Introduction 

CDM models predict two orders of magnitude more dark matter halos than satellites 
observed around the Milky Way (Moore et al. 1999; Klypin et al. 1999). Solutions to this 
problem include self-interacting dark matter (Spergel & Steinhardt 2000), truncated power 
spectra (Kamionkowski & Liddle 2000) and the restriction of gas accretion to the lowest mass 
dark matter halos to be before the epoch of reionization (Bullock, Kravtsov & Weinberg 
2001). The first two solutions would get rid of the smallest dark matter halos entirely, while 
the latter would predict that only one percent of the sateUite halos actually contain stars. 

In this paper, rather than attempting to solve the problem, we discuss how we might 
tell observationally whether it is a problem in the first place. The sateUite dark matter halos 
in the CDM models comprise about 10% of the total mass of the parent galaxy on roughly 
isotropic orbits distributed throughout the galaxy (Font & Navarro 2001). One idea is to 
look for signatures of these lumps around external galaxies in gravitationally lensed images 
of background quasars (Chiba 2001; Metcalf 2001). We concentrate our study rather closer 
to home. If such lumps truly exist around the Milky Way then we would expect them to 
have some dynamical influence on the rest visible Galaxy. For example, according to Lacey 
& Ostriker (1985) if the halo is entirely made of lO^M© black holes they would significantly 
heat our Galactic disk. Font & Navarro (2001) used numerical simulations to ask whether 
the coldness of our disk could be used in a similar fashion to limit the distribution of dark 
matter lumps seen in cosmological simulations that could be orbiting the Galaxy. Using one 
realization of lump masses and orbits taken from a ACDM model of a Galaxy-sized halo, 
they found that the heating caused by the lumps was less than the heating observed in the 
stellar populations in the disk and concluded that the disk is not an efficient probe of this 
model. 

Streams of debris from the destruction of Galactic satellites are another example of 
cold structures within the Milky Way that could be scattered by substructure in the poten- 
tial. These streams tend to align along a single orbit (Johnston, Hernquist & Bolte 1996; 
Johnston 1998; Helmi & White 1999) and hence would individually have lower cross section 
to interactions than the Galactic disk. Nevertheless we chose to examine the response of 
these to perturbations rather than the disk for several reasons: they will explore the outer 
regions of the Galactic potential that the disk does not experience; we expect other sources 
of heating in this region (structures such as the bar, spiral arms and giant molecular clouds 
in the disk) to be negligible; if streams from several satellites could be studied, we have the 
potential of probing a larger portion of the Galaxy; finally, these streams are often even 
colder than the disk itself and hence should be more sensitive to scattering. 

We approach this problem of scattering numerically rather than analytically both be- 
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cause of the nature of the predicted CDM mass distribution of lumps and the low cross 
sections of the streams. Scattering by the few most massive lumps in the distribution is 
expected to be most important, and the exact degree of scattering will depend on the rel- 
ative orientation of the lump and debris orbits. Hence the process is dominated by a few 
strong encounters rather than one that can be modeled (say) by integrating over many weak 
encounters in the impulsive regime. Nor would a simple analytic representation include how 
the wake in the halo, excited by the most massive lumps, would affect the stream orbits. 
Weinberg's studies (Weinberg 1998a, 1995) of the affect of the Large Magellanic Cloud on 
the disk of the Milky Way suggest that the inclusion of such a wake in the halo introduces 
a large enhancement to the expected response of the disk. 

In §2 we present our numerical approach to implementing the experiments described 
above. In §3 we compare the evolution of tidal debris in smooth and lumpy potentials and 
propose an algorithm for distinguishing between the two with currently available data. In 
§4 we apply our results to the stream of carbon stars known to be associated with Sgr. We 
summarize our conclusions and outline future prospects in §5. 

Note that during the final stages of preparation of this manuscript two other papers 
have appeared on this subject (Mayer at al. 2001; Ibata et al. 2001a). We discuss their 
relation to our own work in §5. 

2. Methods 
2.1. Particle Distributions 

We model our system with three sets of particles: 

-^haio smooth halo particles. These are equal mass, and initially distributed as an equilibrium 
Hernquist (1990) model. This model is assumed to represent a 2.3 x lO^^M© halo, with 
scale length Ohaio = 30.6kpc, which mimics the NFW profile (Navarro, Frenk & White 
1996, 1997) of a 1^200 — 200km/s halo at 2; = in a ACDM Universe (taken from Navarro 
& Steinmetz (2000) — paper since withdrawn, see below). Navarro & Steinmetz (2000) 
showed that such halos are incompatible with the apparent mass distribution in the 
Milky Way, but we will nevertheless adopt them since we wanted to place our model 
in a consistent ACDM universe. 

N\ump lumps. These are each represented by rigid Hernquist (1990) models, with initial 
positions and velocities chosen at random from the same Hernquist (1990) distri- 
bution as above. The masses are chosen at random from a power-law distribution 
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dN/ dmiump oc mjyjj^p , whose form is taken from the end point of N-body, cosmologi- 
cal simulations of CDM Universes (Klypin et al. 1999; Moore et al. 1999). The upper 
and lower limits of the distribution are chosen so that 10% of the parent halo mass 
is contained in lumps with masses in the range 7 x 10^ — 2 x lO^^M© which approxi- 
mately matches the level of substructure observed in N-body simulations (Moore et al. 
1999). Note that while figures 4 and 5 in Klypin et al. (1999) suggest that the degree 
of substructure in their ACDM model is similar to their CDM case, Font & Navarro 
(2001) found the numbers of sateUites in their ACDM model to be roughly half those 
seen in their CDM model. We therefore adopt the Klypin et al. (1999) and Moore et 
al. (1999) level of substructure as an upper limit to the degree of heating expected in 
a ACDM universe. 

The scale lengths oiump are set to aiump = Sr^ and the masses miump are chosen so that 
the potentials mimic "Universal" Navarro, Prenk & White (1996) profiles to within 10% 
out to radii lOr^, with the mass-dependence on concentration for a ACDM Universe 
taken from Navarro & Steinmctz (2000). Note that since completion of the numerical 
portion of this project the results of Navarro & Steinmetz (2000) have been withdrawn 
and the models we adopted have been shown to be too concentrated. This will again 
cause a shght overestimate of the efficiency of scattering. 

Hernquist (1990) profiles were used rather than NFW profiles because of the saving 
in cpu-time. We have not tidally truncated the models but do not expect this to 
significantly affect our results. 

Figure 1 summarizes the mass and scale distribution used in the simulations with the 
solid lines representing the parameters for the Hernquist model and the dotted lines 
representing the parameters for the equivalent NFW model. 

Attest test particles. These are massless and arc cither distributed on circular orbits at 
r = 0.5, 1.0, 2.0 and 4.0ahaio) oi' on a range of orbits to mimic debris streams from 
Sgr. The initial conditions for the latter were generated by running a simulations of 
satellites of mass 10^ disrupting in a rigid representation of the halo potential described 
above (for simulation technique see Johnston, Hernquist & Bolte (1996)) along an orbit 
similar to that expected for Sgr (see Ibata et al. (2001a)). The time at which particles 
were lost from the satellite were noted, and the positions and velocities of a set of 500 
particles used as initial conditions, starting from a point shortly after the pericentric 
passage at which they were first unbound. 
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2.2. Force Computations 

The influence of tlie smootli lialo particles on each other and all other particles is 
calculated using a code that employs basis function expansions to represent the potential 
(Hernquist & Ostriker 1992) (we used the MPI version of the code, adapted from the original 
by Steinn Sigurdsson and Bob Leary). This effectively smooths over all strong encounters 
(which we expect to be unimportant for this component), but does follow both large-scale 
collective fluctuations and the halo's response to disturbances by lumps. 

The influence of the lumps on each other and all other particles is calculated directly 
using the analytic form of the Hernquist (1990) potential. 

The test particles respond to both the smooth halo and the lumps, but do not otherwise 
interact. This simplification is valid for our application of tidal streams, but would be a 
much poorer approximation were we trying to represent a system such as a disk which has 
significant self-gravity. 

2.3. Integration 

Simple leap-frog integration was used throughout. The major computational expense 
arose from the large number (A^haio > 10^) of smooth halo particles required (see §2.5). 
However, since we were not interested in modeling strong encounters in this component in 
detail, we were able to use a large timestep {dt — 0.08 in simulation units or about 4 Myears) 
for these particles corresponding to 1/lOOth of the dynamical time at the half mass radius 
of the parent galaxy. The lump and test particles were integrated with a much smaller 
timestep dtm — dt/2"'^'\ where nm was chosen so the increase in computational cost for these 
integrations was not significant. For A^haio = 10^, Ni^^ap = 20 and Attest = 1000, we could 
take niii = 8 (or 256 small steps for every large step and dtm ~ 1.6 x 10'^ years) with just 
a 20% increase in cputime, and we double checked our simulations by rerunning them with 
nm = 9 to confirm convergence of the results. In general, we were looking to accurately follow 
encounters between features of a kpc with lumps traveling at 200km/s so we would expect 
dt\ii — (lkpc/200km/s)/100 ~ 50000 years to be a sufficiently small integration time-step. 

2.4. Initialization 

The smooth halo distributions were allowed to run for ten dynamical times to erase any 
effects from generating the initial conditions. For each set of lumps, the combined smooth 
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halo and lumps were run for an additional ten dynamical times as the lumps were grown 
slowly from zero to full strength to ensure no effects from suddenly introducing them. Finally, 
the output of the latter simulations were used as initial conditions in which the test particles 
were also run. 



2.5. Required Particle Number, Number of Simulations and CPU-cost 

The finite number of halo particles will introduce potential fluctuations in addition to 
those due to the lumps (Weinberg 1993, 1998b) so wc need to ensure that the dominant 
cause of test-particle scattering is due to the lumps and not to numerical noise. To quantify 
the level of fluctuations in the halo due to the finite number of particles wc first ran models 
with logA^haio = 4,5,6 and 7 in isolation for ten dynamical times and recorded the basis- 
function expansion coefficients Anim calculated by the code during this time (where n refers 
to the radial basis function and (/m) refers to the spherical harmonic — see Hernquist & 
Ostriker (1992)). The total potential energy can be written in terms of this expansion as 
$ = Tinim-^liim- Hcuce the potential energy associated with each of the spherical harmonics is 
(^im = TsnA^i^. The points in Figure 2 show the dispersion in the lowest order recorded 
during the isolated simulations. The dotted line shows the same measurement made in a 
simulation containing A^haio = 10'' particles and just one NFW lump. The plot suggests that 
to correctly model the infiuence of lumps on tidal streamers requires A^haio ~ 10'' particles 
since otherwise scattering due to the finite resolution of the halo will compete with the effect 
being measured. 

A second numerical consideration is that, so long as the most massive particles have the 
greatest infiuence, the size of the effect we are trying to measure will depend on the exact 
choice of lump and streamer orbits. To gain a fair impression of the general evolution it is 
necessary to run many simulations with different lump orbits. 

The combined requirements of large particle number and multiple runs supports our 
choice of running idealized simulations of fully formed galaxies over a restricted amount of 
time rather than trying to perform this experiment fully self-consistently in a cosmological 
context, which would add additional computational expense. Each simulation was first 
"relaxed" for 500 large steps and then followed for 1000 large time steps, or ten dynamical 
times (corresponding to about 4 Gyears) for a total cpu-cost of ~ 225(A^haio/10^)hours. We 
ran more than 50 such simulations (a subset of which are presented here) . To reduce this to 
manageable proportions (in memory and time), simulations were run on multiple nodes of 
the Wesleyan Beowulf-class supercomputer (WesWulf). 
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2.6. Parameter Sets 

We present the results from a total of 37 simulations in this paper. Each simulation 
contained A^haio = 10^ halo particles and Attest = 4000 test particles. The "control" simulation 
contained no halo lumps. The remaining 36 consisted of 6 sets of 6 random realizations of 
lump orbits for A^iump = 1, 4, 16, 64, 128 and 256. Each set of simulations contained the same 
lump mass distribution, with increasing numbers of lumps exploring further down the mass 
function. Since test particles were distributed on orthogonal orbits our final results contained 
12 different realizations of debris distributions, all starting from the same initial conditions 
but experiencing different time dependent potentials. 

3. Results 

3.1. Evolution of Orbits in Lumpy Potentials 

As an example of evolution in a lumpy potential Figure 3 plots the phase-space distri- 
bution of particles initially on circular orbits at r = 0.5 and 1.0 kpc at times t = 0.8,2.4 
and 4Gyrs from one simulation containing A^iump = 256 lumps. The left hand panels show 
positions in the orbital plane while the right hand panels show what observations we might 
make of the stars (assuming we could identify the initial orbital plane of the debris): angular 
distance d9 from the orbital plane and line-of-sight velocities v as a function of angle ^ along 
the debris. The two things immediately apparent are: (i) the orbital plane precesses over 
time (visible as a sinusoidal shape in the lowest d9 vs ^ panel); and (ii) particles initially 
distributed smoothly along the orbit become bunched in angular position and velocity. 

Figures 4 and 5 repeat the above plots for the tidal debris particles in simulations with 
zero and A^iump = 256 lumps. Comparison of the right hand panels of these plots clearly show 
the orbital precession described above in the lumpy case. There is also some indication of 
additional structure along the orbit, but this is somewhat hidden by the intrinsic non-uniform 
distribution of debris along the orbit. 

3.2. Interpreting Observations 

From our qualitative assessment of Figures 3-5 we know we need to design a statistic 
that is sensitive to scattering in both angle and velocity along a debris stream rather than 
to large scale effects such as overall heating. (Although precession is an additional affect it 
is not observable since do not in general know the original orbital plane of the satellite.) A 
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natural choice is to look at fourier series in the observed quantities: 

= iSAexp™"'^/''— |. (1) 

In general we expect the low order to contain signal that could be due to periodic nature 
of debris orbits (in v) or debris density (in ^) along the orbit. However, higher Bj^ should 
be sensitive to scattering in the debris distribution. 

To mimic how this might be apphed to a real data set (such as Sgr) we: 

(i) "Observe" the Galactic latitude and longitude and line-of-sight velocity v oi n debris 

particles in our simulations from a viewpoint 8kpc from the Galactic center; 

(ii) Define the orbital plane to be the best fit great circle to the angular data (Johnston, 
Hernquist & Bolte 1996); 

(iii) Find and d6 relative to this great circle; 

(iv) Define a "scattering index" by summing over the fourier expansion - 

B = VSm=5,10S^ (2) 

where 

Note that the limits in the summation were specifically chosen so that the statistic is not 
sensitive to large scale effects such as intrinsic non-uniformity of tidal debris. 

The top panel of Figure 6 plots B^ against m calculated from all 500 of the debris 
particles from the IO^Mq satellite at the end of the simulations which each contained 256 
lumps (open squares and solid lines) and contrasts these with the simulation containing no 
lumps (closed squares and bold lines). The bottom panel plots the scattering index for all 
the simulations as a function of number of lumps (open squares) with the bold line showing 
B for the no-lumps case. In general, we find we expect to be able to use this index to 
distinguish between sAGDM and non-lumpy Milky Way halos roughly ninety percent of the 
time. There is a tendency for the majority of the scattering to be caused by the very largest 
lumps. However, the large variance in the results for a given A^iump emphasizes that the 
degree of scattering is very sensitive to the relative orientation of lump and debris orbits. 

4. Application to Observations: Sagittarius' Carbon Star Stream 

Ibata et al. (2001a) report the discovery of a set of carbon stars which closely align with 
the great circle defined by Sgr's position and proper motion. They find that this stream 
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is too thin to be easily explained in models of the Galaxy that have oblate halos because 
debris orbits in non-spherical potentials process beyond the width of the stream within the 
typical lifetimes of the Carbon stars (47 of the 104 stars lie within 10 degrees of Sgr's great 
circle which intercepts less than one quarter of the total survey area). Such a cold stream 
provides an obvious testing ground for models that predict dark matter substructure around 
the Milky Way 

Carbon stars have ages ranging from a few to 6 or 7 Gyears. Hence we expect the length 
of our simulations (4 Gyears) to be a fair representation of how long this debris may have 
been orbiting independently 

Figure 7 repeats Figure 6 but for 47 tidal debris particles chosen at random at the end of 
the simulations. In this case the particles were again viewed at point 8kpc from the Galactic 
center chosen slightly out of the orbital plane of the satellite to reflect the Sun's orientation 
relative to Sgr's own orbit. To mimic the real survey, we also restricted our simulated survey 
to only those particles further than 30 degrees from the Galactic disk and within 10 degrees 
of the best-fit great circle of the full debris distribution. In the upper panels of Figure 7 the 
points and lines represent the results from just one set of particles in each simulation that 
contained A^iump=256. In the lower panels the points represent the results for one set of 47 
particles chosen from each simulation. The lower shaded region in the lower panels shows 
dispersion around the average of the results for ten different particles sets in the no-lumps 
case. The bold dashed line in both panels is the result of the identical analysis performed 
on the Sgr carbon star data set. The figure suggests that the Sgr carbon star set exhibits a 
level of substructure inconsistent with debris orbits on a smooth, spherical potential. 

The "lumps" that we are already aware of in the Milky Way's halo consist of the Large 
and Small Magellanic Clouds and the eight other dwarf spheroidal satellites. Since our 
results clearly depend on the exact orientation of the orbits of the largest lumps masses 
relative to the debris orbits we ran one final simulations in which we integrated the Sgr 
debris in a halo containing a satellite of mass lO^^M© in an orbit like that of the LMG 
(roughly perpendicular to Sgr's orbit with a pericenter comparable to Sgr's apocenter). The 
results of the application of the scattering index to to the output from the simulation are 
shown as bold dotted lines in Figures 6 and 7. We conclude that for the current Sgr data 
set the degree of scattering is entirely consistent with debris perturbed by the LMC alone. 
Moreover, Figure 6 implies that even with a larger data set associated with Sgr it will be 
difficult to distinguish between a sACDM halo and a smooth one. 
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5. Conclusions and Future Prospects 

In this paper we showed that it is possible to distinguish between smooth and lumpy 
Milky Way halos by quantifying the coldness of tidal streams. We proposed a "scattering 
index" , based on position and radial velocity measurements of stars, that is sensitive to small 
scale perturbations in the debris rather than large scale effects such as variations in debris 
density or hne of sight velocity (due to its basic dynamical properties) . We found that this 
statistic, when applied to measurements of 500 stars in a single debris trail, could distinguish 
between smooth, spherical Milky Way models and those containing a level of substructure 
consistent with sACDM models 90% of the time. Most scattering was due to the few largest 
lumps and the degree of scattering was very sensitive to the exact orientation of lump and 
debris orbits. 

These results agree qualitatively with those of Mayer at al. (2001) and Ibata et al. 
(2001b). Mayer at al. (2001) present images from simulations of tidal tails evolved in a 
fully self-consistent sACDM models to illustrate to what extent the tails are disrupted. 
This approach has the advantage of being able to follow the evolution of debris within a 
cosmological context, but (as noted in §2.5), is limited by the resolution and cost of such 
simulations to single realizations of galaxies that may not be adequately resolved to conquer 
intrinsic scattering due to numerical noise. Ibata et al. (2001a) simplify our own approach 
by modeling only the influence of lumps on the debris, not accounting for the halo wake in 
their simulations and hence (as they note) underestimating the true magnitude of the lump 
influence. Nevertheless, they demonstrate that scattering due to dark matter lumps should 
be easily detectable in the angular momentum distribution of debris from globular clusters 
— an experiment that will become feasible with the launch of the GAIA satellite in the 
second decade of the twenty first century. 

In contrast to Ibata et al. (2001b), our own aim was to look at the feasibility of con- 
straining the halo dark matter distribution with data available today or in the near future. 
Following this aim we applied our scattering index to the carbon star stream associated with 
Sgr and found that it contains more scattering than would be expected if it were orbiting in 
a smooth Milky Way. We then demonstrated that the level of scattering in Sgr's debris is 
entirely consistent with perturbations by the LMC alone. We conclude from our study that 
the current data is unable to place limits on the level of dark matter substructure in the 
Milky Way's halo. Moreover, the specific ahgnments of the LMC's and Sgr's orbits means 
that a larger sample of Sgr's debris is unlikely to improve the sensitivity of our statistic. 

Two other factors could be responsible for the large B measured for the Sgr sample: 
we have not taken into account the contamination of the carbon star sample by non-Sgr 
stars; and we have assumed that the Galaxy is perfectly spherical and that the stream is not 
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broadened by precession of the debris orbits. Despite these hmitations we remain optimistic 
about the use of tidal debris to constrain dark matter substructure in the future. The former 
problem should be solved once we have a large enough data set to identify Sgr debris by 
continuous velocity variations across the sky, and the latter could be addressed with a clean 
sample by looking for anomalous local scatterings in velocity. 

In addition, there is mounting evidence to suggest that many others of the Milky 
Way's satellites (and globular clusters) have associated debris (Grillmair et al. 1995; Ir- 
win & Hatzidimitriou 1995; Majewski ct al. 2000b). This evidence is currently limited to 
ovcrdcnsitics of material close to these objects and hence too young dynamically to be of 
use for our study. However several deep-halo surveys are currently underway that have the 
potential of tracing these streams out further (Ivezic et al. 2000; Yanny et al. 2000; Majewski 
et al. 2000a; Morrison et al. 2000). If extensive debris is found then these satellites, all of 
which have lower masses than Sgr, would provide even colder streams which should be even 
more sensitive probes of substructure. They are also on different orbits, exploring other 
regions of the halo, where the LMC's influence will be less dominant. Looking even further 
to the future, astrometric satellite programs such as NASA's SIM and ESA's GAIA will 
provide the two additional dimensions of proper motion to allow a more rigorous measure of 
scattering to be defined, such as that proposed by Ibata et al. (2001b). 

We would like to thank Mike Irwin for generously making the Sgr carbon star data 
set available to us, Steinn Sigurdsson for sending us the MPI version of the SCF code and 
Ben Moore and Julio Navarro for helpful comments on the sACDM model used. KVJ was 
supported by NASA LTSA grant NAG5-9064 and CH was supported by funds from Wesleyan 
University. DNS was partially supported by NASA ATP grant NAG5-7154. 
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Fig. 1. — Masses, scales and cumulative mass of the top A^iump chosen at random from 
our ACDM spectrum. Sohd lines represent parameters for Hernquist models chosen to be 
equivalent to the NFW models (dotted lines) 



- 15 - 



O 



-1 ^ 

-2 % 
-3 
-4 



-5 
-3 

-4 

-5 

-6 

-7 
-3 

-4 

-5 

-6 

-7 
-4 

-5 

-6 

-7 t 

-8 
-4 

-5 

-6 

-7 

-8 
-5 

-6 

-7 

-8 

-9 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
(l,m) = (0,0)- 



(l,m) = (l,0)- 



(l,m) = (l,l)1 



(l,m) = (2,0) E 



(l,m) = (2,l)l 



(l,m) = (2,2) 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



3 4 5 6 7 8 
log(N) 

Fig. 2. — Dispersion (calculated over 10 dynamical times) in the potential energy associated 
with each {l,m) spherical harmonic in a halo realized by iVhaio particles. The dashed line 
shows the same calculation for a halo with = 10^ particles with the most massive CDM 
lump orbiting in it. 
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Fig. 3. — Final positions (left hand panels) and "observations" (right hand panels) of test 
particles initially on circular orbits after 1.3 (top panels), 2.6 (middle panels) and 4 (bottom 
panels)Gyears of evolution in a halo containing N — 10^ particles and A^iump=256 lumps. 
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Fig. 4. — Final positions (left hand panels) and "observations" (right hand panels) of test 
particles initially on Sgr debris orbits after 1.3 (top panels), 2.6 (middle panels) and 4 
(bottom panels) Gyears of evolution in a halo containing N — 10^ particles and Mump^O 
lumps. 
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Fig. 5. — Final positions (left hand panels) and "observations" (right hand panels) of test 
particles initially on Sgr debris orbits after 1.3 (top panels), 2.6 (middle panels) and 4 
(bottom panels) Gyears of evolution in a halo containing N — 10^ particles and A^iump=256 
lumps. 
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Fig. 6. — Upper panel shows components of the scattering index B^, calculated from the 
final observations of 500 Sgr debris particles, plotted as a function of fourier number m (see 
equation [3] for definition). Open squares and solid lines show results for all realizations 
of halos with Ni^mp = 256 lumps on random orbits. Filled squares and bold hues show the 
result for evolution in a smooth halo. Lower panel shows the scattering index B (see equation 
[2]) calculated for all 500 debris particles at the end of each simulation, as a function of the 
number lumps in the simulation. The bold line indicates B for a simulation containing no 
lumps and the dotted hne is for one containing a single lump of mass lO^^M© on an LMC-like 
orbit. 



-20- 



0.04 




5 10 15 20 



m 



T 




J I I I I I I I lJ I I I I I I I lJ I L 

1 10 100 

log(N) 



Fig. 7. — As Figure 6 but for 47 particles selected to be within 10 degrees of the best fit 
Great Circle to the debris particles at the end of the simulation, and at Galactic latitudes 
l&l > 30 degrees. The bold dashed hues in both panels shows the result of applying the 
same statistic to the 47 Carbon stars found within 10 degrees of Sgr's orbital plane. In the 
lower panel, the lower shaded region indicates the dispersion around the mean for 10 sets of 
47 randomly selected particles from the simulation containing no lumps. The upper shaded 
region shows the same region for a simulation containing just one lump, mass 10^^ Mq on an 
LMC-hke orbit. 



